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ABSTRACT 

We report on high contrast mid-infrared observations of Fomalhaut obtained with 
the Keck Interferometer Nuller (KIN) showing a small resolved excess over the level 
expected from the stellar photosphere. The measured null excess has a mean value 
of 0.35% ± 0.10% between 8 and 11 fim and increases from 8 to 13 microns. Given 
the small field of view of the instrument, the source of this marginal excess must be 
contained within 2AU of Fomalhaut. This result is reminiscent of previous VLTI K- 
band (~ 2/im) observations, which implied the presence of a ~ 0.88% excess, and argued 
that thermal emission from hot dusty grains located within 6 AU from Fomalhaut was the 
most plausible explanation. Using a parametric 2D radiative transfer code and a Bayesian 
analysis, we examine different dust disk structures to reproduce both the near and mid- 
infrared data simultaneously. While not a definitive explanation of the hot excess of 
Fomalhaut, our model suggests that the most likely inner few AU disk geometry consists 
of a two-component structure, with two different and spatially distinct grain populations. 
The 2 to 11 microns data are consistent with an inner hot ring of very small (~ 10 to 
300 nm) carbon-rich grains concentrating around 0.1 AU. The second dust population - 
inferred from the KIN data at longer mid infrared wavelengths - consists of larger grains 
(size of a few microns to a few tens of microns) located further out in a colder region where 
regular astronomical silicates could survive, with an inner edge around 0.4 AU to 1 AU. 
From a dynamical point of view, the presence of the inner concentration of sub-micron 
sized grains is surprising, as such grains should be expelled from the inner planetary 
system by radiation pressure within only a few years. This could either point to some 
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inordinate replenishment rates (e.g. many grazing comets coming from an outer reservoir) 
or to the existence of some braking mechanism preventing the grains from moving out. 

Subject headings: instrumentation: interferometers — infrared: stars — stars: circumstellar 
matter — stars: individual (Fomalhaut) 



Introduction 



Fomalhaut is a bright (V = 1.2) relatively young (~ 200 Myr, |Pi Folco et al.l (120041 )) A4 main 
sequence star l ocated 7.7 pc away featuring a large far infrared (IR) excess first detected by IRAS 
( lAumannlll985l ) . The su spected surrounding debris disk was first resolved i n the far-IR by the K uiper 
Airborne Observatory (IHarvey et al.lll996l ). the n in the sub-millim eter ([Holland et al.lll998l ). and 
more recently in scattered light using HST/ACS (IKalas et al.l 120051 ). The latter optical observations 
show a belt of cold dust concentrated around 140 AU from the star, with an asymmetric structure 



consistent with gravitational sculpting by an orbiting planet ( Kalas et al. 



ALMA 350 GHz observations of the Fomalhaut debris ring (IBoley et al. 



2005 



Quillen!l2006f ). Recent 



20121 ) demonstrate that the 



disk parent body population is 13 to 19 AU wide with sharp inner and outer boundaries, and also 
suggest that debris confined by shepherd planets is the most likely origin for the ring's observed 
morphology. Interestingly, the detection of a point source located at 118 AU from the star - just 
inside the inner edge of the ring, matchi ng its predicted lo cation - was reported in 2008 using HST's 
Advanced Camera for Surveys (ACS) (IKalas et al.l 120081 ) . It was interpreted as a direct image of 
the disk-perturbing planet, and named Fomalhaut b. However this detection was made at 600-800 



nm, and no correspo nding; signatures have been found so far in the near-IR range (IKalas et al.l 12008 



Marengo et al.l 120091 ) , where the bulk e mission of such a planet should be expected, nor at 4.5 jttm 
with Spitzer / IRAC (jJanson et al.ll2012l ). The true nature of this HST detected source hence remains 
unclear at this point, but all current models involve dust. It could either be in the form of a large 
circumplanetary disk around a massive planet, or created via a recent collision b etween two Kuiper 



Belt-like objects of radii about 50 km ( ICurrie et al. 1120121 : iGalicher et al. II2012I ). 



Far inside this ring, within a few tens of AU s from the star, a warm dust component has also been 
detected by Spitzer (jStapelfeldt et al.l 120041 ): the InfraRed Spectrograph (IRS) measured an excess 
continuously increasing with wavelength between 17.5 and 34 microns, while direct images obtained 
with the MIPS instrument (Multiband Imaging Photometer for Spitzer) resolved the circumstellar 
disk down to a wavelength of 24 microns. Both of these observations point to a region of compact 
residual excess emission extending inward of ~ 20 AU, but whose spatial and physical structure can 
not be uniquely determined given the limited resolution of Spitzer. The presence of warm dust close 
to the star is independently confirmed by recent 70 /im i mages of Fomalha ut obtained with Herschel/ 
PACS (Photodetector Array camera and Spectrometer) (lAcke et al.ll2012l ). These spectacular images 
show both the cold outer ring and an unresolved excess source co-located with the central star (i.e 
within 5.7", given Herschel beam size) carrying 50% ± 10% of its flux. 

High accuracy long baseline interferometric observations can provide the resolution and contrast 
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required to probe this central region in greater detail. They can focus on the very inner part of Fo- 
malhaut's environment (within a few AU), probing very different astrophysical scales and conditions 
than those studied by single optical/ IR telescopes . near-IR interfe rometric o bservati ons obtained 



at the VLTI were already reported in a first paper (lAbsil et al.ll2009l . hereafter iPaper J ), concluding 



that a K-band excess of 0.88 ± 0.12% is present, and arguing that thermal emission from hot dust 
grains located within 6 AU of Fomalhaut is the most plausible explanation for the detected excess. 
We present here new data obtained at N band (8 to 13 /im) with the Keck Interferometer Nuller 
(KIN) in August 2007 as part of KIN's commissioning "shared risk" science operations, and in July 
2008 as part of the KIN key science program (PI: Phil Hinz). These near-IR (VLTI) and mid-IR 
(Keck) interferometric measurements are then used in conjunction with spectro-photometric data 
to constrain the physical parameters of the inner few AU of Fomalhaut's debris disk. Finally, we 
discuss a few possible astrophysical scenarios which are compatible with the observed disk and dust 
characteristics. 



2. Nulling Set-Up and Observable Quantities 



The o verall KIN system architec t ure and performance are presented in full detail in recent publica- 
tions (jColavita et al.ll2008l l2009t ISerabyn et al.l 120121 ). In brief, four beams are recombined by the 
KIN system. A split mirror located downstream of each Keck telescope adaptive optics system — 
close to a pupil plane — divides the light gathered by each telescope into "left" and "right" beams. 
Interferometric nulling occurs separately between the two Keck left beams, and between the two 
right beams. The resulting nulled output fields are then coherently recombined using a standard 
Michelson interferometer, called the "cross- combiner" . As the optical delay is rapidly scanned inside 
the cross-combiner, one first measures the cross combiner fringe amplitude at null in each of ten 
independent spectral channels covering the full N band (8 to 13 /im), and then "at peak". The null 
depth is defined as the ratio of the cross-combiner fringe amplitudes obtained at null and at peak. 
The rationale for this complex 4-beam combination and modulation, is that the resulting measured 
null depth is free of slow drifts in the incoherent background, a source of strong potential bias for 
ground-based interferometric observations in the thermal IR. Two different scales and baselines are 
then involved: the interferometric nulling baseline of length B ~ 85 m, separating the telescopes 
centers, and the short cross combiner baseline b ~ 4 m, characteristic of the interference between the 
"left" and "right" parts of a given Keck telescope. 

For a perfect instrument, defined as providing a null depth of zero on a point source, one can relate 
the measured monochromatic astrophysical null A ast to the sou rce brightness distrib ution on the sky 
1(6). The observed null can be expressed in the following way (ISerabyn et al.ll2012l ) 



iVast(A) 



/ 1(9) sin 2 (7r£ • 6/X)) cos(2tt6 • 9 / 'X)\l 1 T L (9)T R (9)d9 



(1) 



/ 1(9) cos 2 (nB ■ e/X) cos(2tt6 ■ 9 / 'A))y ' T L (9)T R (9)d9 
where T^(9) and T R (9) designate the sky transmission patterns of the left and right Keck beams, 
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respectively. They are computed from the telescopes orientation and from the overall beam train 
propagation, which includes an intermediary focal plane pinhole. As a result, the field of view FWHM 
is at maximum 450 mas, along the direction perpendicular to the left-right split, corresponding to 
about ± 2AU at Fomalhaut's distance. As shown in Equation [TJ in the case of an extended source, 
the measured null level is not only affected by the long baseline nulling pattern (fast oscillating 
squared sine term), but also by the cross fringe pattern (slowly oscillating cosine term) and by the 
lobe antenna of each single beam. As representative examples of these 3 contributors to the effective 
KIN sky transmission, the left and right panels of Fig. [TJ show the monochromatic (10 /jm) KIN's 
transmission patterns corresponding to the 2007 and 2008 observations of Fomalhaut, respectively. 
Both are derived at the time of meridian transit and mostly differ in the orientation of the short 
baseline (see section 3). Also indicated in these figures is the or ientation of Fomal haut's outer debris 
disk major axis (156° East of North), as imaged by HST/ACS fjKalas et al.ll2005l ). For both epochs, 
the projected Keck to Keck interferometric baseline orientation is ~ 73 degrees away from the outer 
dust disk main axis. 

In the case where no extended emission is present (source angular size <C A/6), the astrophysical null 
expression (Equation [T|) can be approximated by: 



AU(A) 



J I(6)sm 2 (7rB ■e/X)d6 



(2) 



In particular, for a naked star represented by a uniform disk (UD) of diameter 6* (<S X/B), the 
observed astrophysical null is given by: 



iV«t(A) 



7T 



BQ„ 



1A 



(3) 



For a more realistic model, in which a naked star is represented by a limb darkene d disk of diameter 



Alp, with a linear limb darkening coefficient u\, the observed astrophysical null is ( lAbsil et al.l 12006 
201 lh : 

7nB0 LD V {■, 7 U \\ m_a n 

3 



Wast(A) 



4A 



15 



(4) 



In practice, the null depth measured on a point source is not zero but equal to the instrumental null, 
noted iVi ns . Consequently, the observed null depth 7V bs is not strictly equal to the astrophysical null 
iV ast characteristic of the source, and one measures instead: 



AUs = (iV ast + AU/U + ^astiVh 



(5) 



Full details on the t erms contributing; to the instrumental null can be found in IColavita et al.l (120091 ) 



and ISerabyn et al.l (120121 ). For all the observations considered here, both iV ins and iV ast are small 



1 A curious effect of the KIN 4-beam combination is that for emission sources extending further out than A/6 in 
the direction of the cross-combiner baseline b (b ~ 4 m), some regions will contribute "negatively", i.e. effectively 
decrease the observed null depth. This is illustrated by the areas of negative transmission shown in Fig. Q] 
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(< 2%), so that one can use the approximation: 

N ohs = + N ins . (6) 

As in classical stellar interferometry, the instrumental null iV ins is derived from nulling observations 
of calibrator stars with known diameters and limb darkening properties, i.e., with predictable astro- 
physical nulls (Equation HJ). 



3. Observations and Data Analysis 



The b asic observing block of the KIN is a 400ms long "null/ peak micro-sequence" (see lColavita et al. 



( 120091 ). figure 3). Each micro-sequence yields an individual null depth estimate defined as the ratio 
of the cross- combiner fringe amplitudes measured at null and at peak. A null measurement sequence 
consists in 1000 consecutive micro- sequences, from which a mean null estimate is derived, together 
with its standard deviation. 

Six such null measurement sequences were recorded on Fomalhaut in 2007: one on August 28 and five 
more on August 30. Only the six spectral channels with wavelengths shorter than 11 /im provided 
adequate signal to noise and are discussed here. Eight more individual null measurement sequences 
were then obtained on July 17 and 18, 2008. This time, the full 8 to 13 fim range was covered using 
ten spectral channels. 



3.1. Calibrators 



Fomalhaut's observations were interleaved with measurements of four nearby calibrator stars, whose 
characteristics are given in Table [TJ Their limb darkened (LD) diameters are estimated from V and 
K magnitudes (c o rrecte d for interstellar reddening) and surface brightness relationships developed 
by iDi Benedetto! (120051 ) . except for the redder source HD 214966. In this case, the V-K color falls 
outside of Di Benedetto' s relationships applica bility range, and we use instead the surface brightness 
relationships derived by iBonneau et al.l (120121 ) . 

At the central wavelength of each spectral channel, the expected calibrators astrophysical nulls 
are then de rived using Equation HJ assuming a constant N-band linear limb darkening coefficient 
un = 0.12 (jTango &: Davisl 120021 ) for all four calibrators. Limb darkening corrections are expected 
to be small for such stars at 10 /im, and fairly constant over the range of T e g and log(g) covered by 
our calibrators. 

While Di Benedetto (2012) reports relative uncertainties as low as 2% in his LD diameter estimations, 
we use larger (1-cr) error bars to account for uncertainties in K band photometry or mid-IR LD 
coefficients. For HD 222547 and HD215167, which have accurate K band p hotometry, we adopt a 
relative error of 6%. For HD 214966, we use the 7% relative error quoted by lBonneau et al.l (120121 ). 
Finally, HD 210066 only has 2MASS (saturated) photometric measurements available at K-band, 
and we adopt a conservative LD diameter error of 10%. 
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3.2. Calibrated Astrophysical Nulls 

The instrumental null at the time of a given Fomalhaut observation is computed by interpolation 
between the instrumental nulls derived on adjacent calibrators. At each wavelength and for each 
projected baseline length, one computes Fomalhaut's calibrated astrophysical null from the observed 
nulls using Equation [6j The resulting calibrated astrophysical nulls are summarized in Table [2j 

Overall, the length of the long baseline B projected onto the sky plane (B p ) varied between 55.12 m 
and 77.81 m during the observations, while the projected baseline azimuth varied only slightly (be- 
tween 4 7? 4 and 49?9). 

The projected length of the short baseline b, on the other hand, does not change: by construction, 
it is always located in the (pupil) plane perpendicular to the line-of-sight (see Fig. [2]). We note that 
the short baseline orientation was changed by 90° between the 2007 and 2008 observations. In 2007, 
it was located in the plane defined by the zenith and the line-of-sight, while in 2008, it was located 
in the local horizontal plane (see Fig. [2]). These variations in the short baseline azimuth do not have 
a large influence on the measurements but were included in the null computations (Equation [T]) and 
modeling. 



3.3. Null excess leakage 

The next step in the data analysis is to evaluate the fraction of Fomalhaut's calibrated astrophysical 
null that actually comes from the circumstellar environment. This requires the computation of 
the null depth expected from the photosphere alone (naked star scenario, Equation HJ). We use 
Fomalhaut's limb darkened diameter #ld = 2.223 ± 0.022 mas deduced from the latest VLTI K band 



measur ements (|Paper 1) and a constant linear limb darkening coefficient un = 0.06 between 8 and 



13 jttm (ITango fc David 120021 ). For the baselines considered here, the stellar contribution to the total 
null depth typically ranges from 0.5% and 0.2% for wavelengths from 8 to 13 /im. The difference 
between Fomalhaut's observed calibrated null and this purely photospheric leakage is noted "null 
excess" hereafter. It characterizes the contribution of any source of circumstellar emission located 
within the nuller field of view. 

For each of the 14 independent calibrated observations reported in Table El we compute the resulting 
null excess as a function of wavelength. Since there is very little variation in azimuth over the 
full data set but some projected baseline length variations and some changes in the instrumental 
set-up between 2007 and 2008, we further regroup the measurements according to mean baseline 
length ('"short" or "long") and year of observation. This allows us to generate 4 final data sets from 
the original 14, which reduces null measurements uncertainties and the modeling computation time 
(section 4). For instance the 2007 "short" baseline data summarize measurements from the three 
shortest baselines, ranging from 58 to 63 m. From that sub-ensemble of measured null excesses and 
associated error bars, we compute the weighted mean null excess and its 1-cr uncertainty (we use 
the weighted standard deviation). This provides an accurate estimate of the null excess at the 2007 
short observing baseline, which has a mean length of 60.72 m and an equivalent azimuth of 49.77°. 
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Fig. [3] (top left) shows Fomalhaut's measured null excess at this baseline, as a function of wavelength. 
Fig. [3] (bottom left) shows similar curves for the 2007 long baseline data (mean length of 71.88 m, for 
an equivalent azimuth of 48.66°). Analogous results are presented for the 2008 observations in Fig. [3] 
(right panels), grouping the four shortest baselines (55 to 62m, mean length of 58.43m, azimuth of 
49.42°, top right inset), or the four longest ones (73 to 78m, mean length of 75.79m, azimuth of 
47.40°, bottom right inset). These figures show that our 2007 and 2008 data sets are consistent with 
each other (generally to within la), although the short wavelengths excess null looks slightly larger 
in 2007 than in 2008. 

Based on the results presented on Fig. [3] and taking all 2007 and 2008 observations into account, 
the (weighted) average null excess measured bet ween 8 and 11 /im is .35%. Small systematic errors 
have been previously identified in the KIN data (IColavita et al.ll2009f ). potentially biasing the nulls 
measured at different wavelengths in a systematic way during a given night. As a conseque nce, we 
use for our measurements the 1-cx uncertainty level recommended in IColavita et al.l ( 120091 ). which 
represents our best understanding of the instrument. For KIN measurements of stars as bright as 
Fomalhaut in the 8 to 11 /im range, this uncertainty corresponds to 0.2% rms per night, or 0.1% over 
the 4 nights of observations reported here. Our final estimate of Fomalhaut's averaged null excess in 
the 8 to 11 /im region is then 0.35% ± 0.10%. Finally, we note that for the 2008 data - covering the 
whole 8 to 13 /im atmospheric N band window-, the observed excess increases with wavelength, for 
both baseline ranges considered. 



4. Modeling and Interpretation 



Although each of our four data subsets is formally consistent with pure photospheric emission to 
within 2-3 a, an excess null is measured at all wavelengths and for all four of the observing nights, 
covering years with slightly different KIN instrumental set-ups. This suggest the need to go beyond a 
simple photospheric model to explain the observation s. Owing to the various obs e rvatio n al evidence 



for wa rm emission in the Fomalhaut inner system (IPaper II . IStapelfeldt et al.l ( 120041 ); lAcke et al. 



( 120121 ) ) . we hereafter model the measured excess null using purely morphological debris disk models. 
In particular, we do not favor the assumption that either the observed KIN or the VLTI /VINCI 
emission is due to gaseous free-free emission from a stellar wind, and there are three main reasons 
for that: 



For A stars, the photospheric emission still goes as v 2 in the near-IR, while the i/ 0,6 spectral 
slope for free-free emission starts to break down in the mid-IR, where free-free emission becom es 



optically thin, and the slope flattens for shorter wavelengths (e.g. iWright fc Barlowl (119751 )). 
As a result, the free-free emission level relative to the star is expected to be smaller in the 
near-IR than in the mid-IR. Free-fr ee emissi on would probably not be sufficient to reproduce 
the observed VINCI K-band excess (IPaper J ). And if it were, then it would be stronger in the 
mid-IR, which is not compatible with the low level of emission seen by the KIN. 



Second, the size of free-free emission regions decrease at shorter wavelengths. Actually, at the 
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wavelength where free-free emission becomes optically thin, the size of the emitting region is 
of the order of the stellar size. Therefore, we would not expect the near-IR free-free emission 
to be significantly extended (while the previous VLTI observations show the K-band emitting 
region to be at least 10 times larger than the photosphere). 



Third , stellar models predict very small mass loss rates for A-stars (e.g. ,10 M & / yr, iBabel et al 



( 119951 )). which are not compatible with strong free- free emission. 



4.1. Modeling the Mid-Infrared Data with a Solar Zodiacal Disk Model 



Since the observed null depth is a function of source brightness distribution and baseline orientation 
(see Equation [T|), we first need to define a morphological model for any excess emission before 
computing its associated photometric flux. In other words, the true astrophysical excess can not 
be uniquely derived from the data, unless some assumptions are made on its spatial brightness 
distribution. As a first rough attempt to estimate the luminosity of the marginal excess detected, 
we have used a scaled model of the solar system's zodiacal disk around Fomalhaut. This model 
is based on the paramet r ic des cription of the zodiacal cloud observations by the COBE/DIRBE 



instrument iKelsall et al.l (1l998h and is implemented in the Zodipic packageo The dust density 
and temperature profiles are assumed to follow t he power laws deri ved in the solar system case: 



n\r) oc r 



-1.34 



and T(r) oc L s J 2 /r s : with 5=0.467 ( IKelsall et al.lll998l ) and Fomalhaut's luminosity 



L = 17.7 L© ( )Di Folco et al.l 120041 ) . The dust is assumed to extend from the sublimation radius 
corresponding to a temperature of 1500 K (typical of silicate grains found in the solar system), all 
the way to 10 AU, i.e. much further out than the KIN field of view. The exo-zodi disk is modeled with 
the s ame inclination (65.6°) and position angle (156°) as Fomalhaut's outer dust disk ( IKalas et al. 



20051 ). The result of the best fit is a luminosity (and density) scaling factor of 350 with respect to the 
solar case, with 1-er uncertainty of about 100. Such a disk would produce a total emission of 0.55 Jy 
across the KIN field-of-view at 10 /mi, i.e., a flux ratio of ~ 3% with respect to the photosphere. 
This illustrates the fact that due to the nuller sky transmission pattern, the true 10 /im astrophysical 
excess around Fomalhaut could be significantly larger than our observed null excess of ~ 0.35%. This 
is especially true in the case of Fomalhaut, where the disk PA is almost perpendicular to the long 
Keck to Keck baseline (Fig. [T]). 

While such a simple zodiacal disk model can reproduce the full KIN dat a set rea sonably well (Fig. [3]), 
this result seems at odds with the significant K-band excess reported in lPaper II . which would require 
an equivalent dust surface density 5000 times larger than the solar zodiacal cloud to be reproduced. 
This is an order of magnitude larger than the level derived to fit the KIN data. This inconsistency 
forces an examination of dust disk morphologies other than the standard zodiacal disk model, which 
fails to explain the near- and mid-IR observations simultaneously, and which has no real physical 



2 Zodipic is an IDL program for synthesizing images of exozodiacal clouds that can be downloaded at 
http : //asd. gsf c .nasa.gov/Marc . Kuchner/home .html 
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grounding in the case of Fomalhaut. Given the small excess null detected by the KIN (or even 
considering its 3-cr upper limit of just 0.65%), it is clear that the grains responsible for the near-IR 
emission are not contributing much at 10 /im. We will now look for models consistent with this 
observational result. 



4.2. Combined Modeling of Near- and Mid-Infrared Data 



We used the GRaTeR code ( lAugereau et al.lll999l ) to compute a grid of models that we compare with 
our full data set consisting of spectro-photometric and interferometric data gathered in the near and 
mid-IR (see Table. E]and Fig. @J. GRaTeR calculates models for optically thin disks. It is designed to 
simulate spectral energy distributions (SEDs), images and interferometric data with parametric grain 
size and radial distributions, or distributions from dynamical simulations. Both the scattered light 
and the continuum emission of dust grains in thermal equilibrium with a star are computed, using 
the Mie t heory and the Bruggeman effective medium method, depending on the material optical 
constants (IBohren fc Huffman! Il983h . Particular care is given to the removal of grains when their 
temperature exceeds the material sublimation temperature. When multi-material grains are used 
with various sublimation temperatures T su b, the first material to sublimate is replaced by porosity 
(affecting the optical properties with respect to compact grains of the same size). We compute a 
grid of models where we let the following parameters vary: (1) the geometry of the exozodiacal 
disk (assuming azimuthal symmetry before sky-plane projection) defined through its surface density 

/ \-l/2 

profile E(r) = S v2 ( (^)~ 2ain + (^) _2a ) where S is the density at the peak position r , a- m 
the inner slope (fixed to +10 to mimic a sharp inner edge), a the outer slope, (2) the dust size 
distribution (parameters k and a m[n of the classical power-law dn(a) oc a K da valid for grains from 
a m j n to a max = 1 mm), (3) the dust composition Vc/vsi, that assumes mixtures of carbonaceous 
material (volume fraction vc) and astronomical silicates (vsi), (4) the total disk mass A^ust in grains 
up to 1 mm in radius. The range of values used for each parameter is summarized in TabJU This 
parameter space leads to a wide range of different models that we compare w ith the data using bot h 
a classical \ 2 minimization and the statistical Bayesian method described in iLebreton et al.l ( 120121 ) . 
All parameters are assumed to have uniform prior probabilities, at the exception of the treatment 
we use to account for the inner sublimation radius r su b (a, material). In that case, we define prior 
probabilities in order to eliminate all models for which ro < min(r su b). 



4-2.1. Fitting data obtained at all wavelengths 

We first attempt to fit all the near and mid-IR data simultaneously (34 data points, see Tab. [3]). 
The shape of the exozodi spectrum (Fig. H] ) requires very hot grains, which is achieved when these 
are small and close to the star. Indeed, the probability curves derived from the fit Bayesian analysis 
and presented in Fig. [5] (red curves) shows that, when considering each parameter independently, 
the most probable models are found for very small grains (a m i n < 0.08 /im and k < —5.3 with 1-er 
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confidence), confined very close to the sublimation distance (0.07 < tq < 0.14 AU, a < —5.0) §. The 
model requires that these high-temperature grains include a large fraction of carbonaceous material 
> 10), to reduce the ten micron emission feature that submicron-sized silicate grains would 
produce, and to reach higher temperatures. Indeed, the subl imation temperature o f carbon is higher 
than that of silicates (T su b(C) = 2000 K, T su b(Si) = 1200 K. Kobayashi et al.ll2009l ) allowing carbon 
grains to survive closer to the star. The smallest sublimation distances are obtained for the largest 
grains, with min(<i sub (C)) = 0.07 AU, and min(<i su b(Si)) = 0.21 AU (see sublimation curves in Fig. [6] 
for details). 

A clear secondary peak in the disk inner edge probability curve reveals that a second family of solu- 
tions is of statistical relevance, one that uses grains located further out from the star (tq ~ 0.4 AU). 
In fact the model with the smallest overall y 2 (see Tab. [5]) is found among this second family of 
solution^], where ~ 0.2 //m grains orbiting in the ~0.4 AU region, dominate the emission, with 
a characteristic temperature of 1600 K. The disk is largely dominated by thermal emission over 
scattered light at the wavelengths considered here. The result of this fit to the null and spectropho- 
tometric data is shown in Fig. [4] (red curves). The SED is well fitted for a total dust mass of 
~ 4 x 10~ 10 M ffi . The flux contributed by grains of different sizes depending on their distance to the 
star is illustrated in Fig. [6] (upper panels). From Fig. 0] (red curves), it is clear that this overall best 
fit model fails at reproducing the rising null depth observed toward the mid-IR (A > llfim) with the 
KIN in 2008. The model predicts instead a monotonically decreasing contrast toward longer wave- 
lengths and overpredicts the shortest wavelengths null depth from 2008, while staying compatible 
with the 2007 data. Overall, it appears impossible to fit all the data with a single- annulus / single 
grain population model, since the VLTI detection requires very hot dust close to the sublimation 
radius, while the apparent rise of the KIN nulls beyond 11 microns is reminiscent of colder dust. The 
binomial distribution of the best fit inner radius (r , Tab. E]) is also suggestive of a two component 
dust distribution. 

In an attempt to explore this two dust component scenario with the GRaTeR code, which is currently 
limited to a single dust population (already exploring a large 5-parameter space, see Tab]4]), we now fit 
separately the data shortward of \\[im (hereafter labelled "SHORT") and the data longward of 11/im 
(hereafter labelled "LONG"). An improved version of the GRaTeR code, including self consistent 
modeling and radiative transfer calculations through multiple dust components is currently under 
development (Lebreton in prep.) but beyond the scope of this paper. 



3 Formally, these probabilities must read as the probability of the models knowing the data. It must be noted that 
for some parameters, the probability curves tend to peak towards the limits of the parameter space, in a physically 
unmcaningful manner. Thus the confidence intervals we provide must be taken cautiously. In particular, n and a ou t 
should be considered as weakly constrained. 

The most probable solution for individual parameters can indeed be distinct from that found in the 6-dimensions 
parameter space. 
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4-2.2. Fitting data subsets 

The shortest wavelengths (A < 11 /im, observations (26 data points, see Tab. G2) can be fitted by a 
family of models comparable to the one previously discussed: the probability curves (Fig. E]) reveal 
that very small grains located very close to the sublimation limit are favored. The solution now 
has ro in the range 0.08-0.11 AU, i.e. located as close as possible to the star, next to the size- and 
composition-dependent sublimation radius. The best fit model is presented in Tab.[5]and Fig. H] (blue 
curves). Fig. (middle panels) reveals that this solution also consists of sub- micron grains with a min 
in the range 0.01-0.21 \im. Thus both the near-IR and A < 11/xm data can be fit by the same hot 
dust ring located at very small radii. 

Now fitting only the longer wavelengths observations (A > 11/im, 11 points including Spitzer/IRS 
upper limits, see Tab. [3]), significantly different results are found. As can be seen in Fig. @] (green 
curves) and Fig. [6] (lower panels), the KIN observations are well reproduced if grains are located in 
a ring located further out from the star, and two families of solutions emerge with r ~ 0.45AU 
or r ~ 1.0 AU. With an unconstrained k, and a min close to the radiation pressure blowout size 
(obiow(C) = 3.5/xm, abiow(Si) = 2.3/im), the size distribution appears compatible with the dynamical 
and collisional constraints. Any carbon to silicates ratio is equally probable. A wide range of disk 
masses is allowed, from 10 -9 M e for the first family of solutions, to 10~ 6 Mq. The density profile 
does not necessarily need to be very steep, a -1.0 power-law index for instance yields a satisfactory 
fit with a mass of 4 x 10~ 7 M ffi . While it fails to reproduce the observed K-band excess satisfactorily 
(Fig. |4]SED green curve), this model clearly provides the best fit to the overall KIN data. 

Simulating the emission of a fully self-consistent two-rings structure based on a combination of the 
"SHORT" and "LONG" wavelength models presented here is beyond the scope of this work. Yet, 
the rising slope of the nulls depths measurements produced by the LONG model gives us confidence 
that a two-peak model is the best solution to reconcile all of the data. 



5. Summary 

The near and mid-IR interferometric data for Fomalhaut are consistent with a hot debris disk residing 
interior to the habitable zone of the Fomalhaut system (which extends from ~ 4 AU to ~ 6 AU given 
the star's luminosity). The description of this hot debris disk as a single population of grains is not 
sufficient to explain both the near-IR excess flux found by VLTI/ VINCI, and the small mid-IR flux 
reported here. A possible explanation is that the exozodi has a more complex geometry, for instance 
a double-peaked structure. 

In order to be consistent with the KIN measurements, the VLTI near-IR excess requires a large 
population of small (~ 10 to 300 nm) hot dust located very close to the rim defined by the sublimation 
distance of carbon. Such grains are however much smaller than the radiation pressure blowout size 
and they should be pl aced on hyperbolic o rbits and ejected very quickly from this region. A similar 



paradox was found by iDefrere et al.l (120111 ) for the disk of Vega, and unveiling the origin of this hot 



dust is still a challenge to debris disk science. The main difficulty is to explain how such high levels 
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of sub-micron grains can be present, when these grains should be blown out by radiation pressure 
on very short timescales, typically of the order of a dynamical timescale, i.e., much less than a 
year in these inner regions. Is it because these grains are prevented from moving out by a braking 
mechanism, and if yes, which one? gas drag? collisions in a very dense and radially optically thick 
disk? Or is it because the production rate of these small particles is so high in these regions that, 
despite their fast removal, a significant amount is always present at a given time? These issues will 
be explored in a forthcoming paper (Lebreton et al., in preparation). 

On the other hand, the KIN observations detected a small mid-IR excess which appears to increase 
with wavelength between 8 and 13 microns, which constrains the dust location to lie in a colder 
region where astronomical silicates could survive. While the VLTI excess only probes the inner 
hotter component, this second dust population is traced by the KIN data. With micrometer sizes 
compatible with a classical-collisional equilibrium, mid-IR emission at the level seen, would require 
an independent dust population, originating from a much more massive population of planetesimals, 
comparable to the Solar System Main Belt. The low statistics associated with this warm component 
do not allow firm conclusions on the parameters of this grain population, but it is a serious hint that 
a secondary zodiacal belt lies within the field of the interferometer, i.e. inside a few AU. 

Finally, we emphasize that our model has some limitations and should not be considered a definitive 
explanation to the hot excess of Fomalhaut. First, because each component contributes to both sub- 
sets of data, a proper parameterization of the double-population should be used to self-consistently 
describe the exozodi. Second, our treatment of dust sublimation is still ve ry coarse as i t elim inates 
instantaneously the grains when they exceed the sublimation temperature. iKama et al.l ( 120091 ) show 
that a proper treatment of time-dependent sublimation physics allows grains hotter than the cur- 
rent T su b to survive a certain time before vanishing, depending on their size. Allowing for these 
transient grains would likely result in putting some large grains inside the current sublimation rim, 
thus impacting the best-fit models. In a future study, we will improve these aspects of our model, 
hopefully reconciling the results with theoretical predictions. Of course, another possibility is that 
the disk has suffered important variability in the time interval between VLTI and Keck observing 
runs. More observations are needed to answer this open question, and more generally constrain hot 
dust transience around mature stars. 
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Fig. 1. — Left panel: effective Keck Nuller sky transmission at 10 /im, when observing Fomalhaut 
at meridian transit on August 30, 2007 (Julian date: 2454342.95, projected baseline length: 67.6m, 
azimuth: 48?7). North is up, East is to the left. High frequency fringes correspond to the long baseline 
separating the telescopes. The low frequency modulation is produced by interference between the 
sub-apertures of a single Keck telescope ("cross combiner" fringes); these fringes are aligned with the 
North-South direction when observing a star at transit. The contours indicate inner regions of the 
Fomalhaut system (i = 66°, PA = 156°), showing that the KIN is sensitive to dust emission in the 
0.05 to 2 AU range. Right panel: Keck Nuller sky transmission at 10 /im, when observing Fomalhaut 
at meridian transit on July 17, 2008 (Julian date: 2454664.07, projected baseline: 67.8 m, azimuth: 
49?6). For the 2008 observations, the telescope pupil was rotated by 90° with respect to 2007. As 
a result, the low frequency ("cross combiner" fringes) are now aligned with the East-West direction 
when observing a star at transit. 
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V 



Fig. 2. — Configuration of the long baseline B (left) and of the short baseline b (right) of the 
Keck Interferoemter Nuller relative to the Fomalhaut planetary system (represented by its cold dust 
annulus, with a position angle of 156°). All the dashed lines and curves are located within the same 
plane, including the zenith and the line-of-sight. The two bold dots represent the centers of the two 
sub-pupils defined on a given Keck pupil. 
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Fig. 3. — Left panel: 2007 calibrated excess null depth measurements of Fomalhaut plotted as a 
function of wavelength (crosses with error bars). Data were obtained at six different baselines and 
grouped into one short equivalent baseline (~ 61m) and one long (~ 72m, see text for details). 
Right panel: 2008 calibrated excess null depth measurements of Fomalhaut plotted as a function of 
wavelength (crosses with error bars). Data were obtained at eight different baselines and grouped 
into one short equivalent baseline (~ 58m) and one long (~ 76m, see text for details). In both 
panels, the expected photospheric null depth has been subtracted from the original KIN data to 
construct the "excess null", which reveals a possible circumstellar excess. The excess null created by 
a 350-zodi exozodiacal disk is shown with triangles and dashed lines for comparison. 
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Fig. 4. — Upper 4 panels: KIN calibrated null measurements (purple diamonds) and results of best 
fit models (Tab. triangles). "ALL" wavelengths model (red curve): uses all interferometric and 
spectrophotometric data from 2 to 13 /im. "SHORT" wavelengths model (blue curve): uses 2 to 11 
[xm data only. "LONG" wavelengths model (green curve): uses 11 to 13 /im data only . Bottom panel: 
solid black line: synthetic photosphere SED model. Purple diamonds: mid-IR spectro-photometric 
measurements and K-band VLTI excess. Red, blue and green curves: resulting emission from the 
disk, for each of the 3 best fit models (same color codes as for upper panels). Solid thick lines: total 
(scattered + thermal light) emission. Dotted lines: thermal light only. A summary of data points 
used for each of the 3 models can be found in Tab. [3j 
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Fig. 5. — Probability of the models knowing the data calculated over a grid of ~ 4, 000, 000 solutions 
for the fit to all observations (ALL), the short wavelengths observations (A < 11/im, SHORT), and 
the long wavelengths observations (A > 11/im, LONG). Depending on the exact composition, the 
sublimation distances range from approximately 0.07 AU to 0.2 AU for carbon and from approxi- 
mately 0.2 AU to 0.8 AU for silicates. It is remarkable that the probability curves for the parameter 
tq is encompassed within these intervals 



1.0 



0.6 



0.2 [ 



0.0 . 



10 



15 




0.01 



0.10 



1.00 10.00 



100.00 




1.0 



t5 0.8 



0.6 



- ■ ALL 
■■■■ SHORT 

- - LONG 



L -J 



-6.0 -5.5 -5.0 -4.5 -4.0 -3.5 -3.0 -2.5 




-20 - 



Fig. 6. — Flux density maps as a function of grain size and distance from the star. Left: A = 2 /xm, 
right: A = 13 fim. At each wavelength, maps are shown for 3 models (Tab. [5], either considering 
ALL wavelengths data (upper panel), SHORT wavelengths data only (A < 11/xm, middle panel), or 
LONG wavelengths data only (A > 11/xm, bottom panel). For a given model, the 2 and 13 fim flux 
densities essentially differ by a scaling factor, except that larger grains contribute slightly more to 
the relative emission at 13 fim than at 2 fim due to their lower temperature. The maps are weighted 
according to the real (fitted) size distributions and they corresponds to the total flux per unit grain 
size per elemental annulus of radius distance: the flux in janskys can thus be retrieved through direct 
integration over distance and grain size. Solid lines: sublimation distances of silicates (green), and 
carbon (when relevant, white). The contours are plotted every power of 10 between 10 _1 and 10~ 7 
of the maximum flux density. 
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Table 1. Calibrators Characteristics 



Star 


Type 


Dec 


R.A. 


6> LD (mas) 


HD 222547 


K4/5III 


23 41 34 


-18 01 37 


2.70±0.16 


HD 210066 


M1III 


22 08 26 


-34 02 38 


2.58±0.26 


HD 214966 


M5III 


22 42 22 


-29 21 40 


3.89±0.27 


HD 215167 


K3III 


22 43 35 


-18 49 49 


2.67±0.16 



Note. - Calibrators used for the Fomalhaut observa- 
tions. 



Table 2. Summary Of Fomalhaut KIN Observations 



Date 


MJD 


B f , 


Az 


iV 8 .25 


A 8 .74 


A9.20 


A9.75 


A?10.22 


^10.68 


iVll.21 


iVli.72 


A r 12.20 


Ni2.69 


no /oo 1 CV7 
Uo/ zo/U f 


5434U.4/ foY 


62.54 


49.88 




i\\ n_i_ i oo 


Q A OJ_ 1 A n 


-i r\n'C_i_ 1 n 1 


1 0_L 1 70 


on f; _lo q o 










no / on /n7 

Uo/ou/u r 


54342. 4U3DO 


76.14 


47.41 




r/in _i_ i a q 
o4U±14o 


1 CY7 A _l_ 1 cn 


i o a 1 _i_ ion 


1 1 Q 7_LOflO 


n o o _i_ o o a 










no 1 on 1 r\ r 7 

Uo/oU/U7 


54342. 42594 


72.46 


48.78 


£?nn I 1 a a 


T1 T 1 1 OC 

717±loo 


o>7/i i no 


C7n I 1 

87y±17o 


ono I 10/ 


i n tr 1 I oo / 










no /on /n , 7 
Uo/oU/U ( 


D4o42. 45259 


67.10 


49.75 


T CT O _J_ "I C A 


coo _i_ it:/? 


cm _i_ i ,i c 
uzi±14o 


a i n_i_ i "70 


rn _l on 1 


O Q _LOC /I 










no /on /n^ 
Uo/oU/U ( 


54342.4/ /8/ 


61.10 


49.81 


i rrn_i_ i no 


,i nnxi oo 


uUU±lli 


O 1 O _L_ "1 O "1 


a to _i_ ion 


rril J_ O O O 










no /on 1 r\ r 7 

08/30/07 


54342.48829 


58.59 


49.54 


coi 11 on 

o31±129 


t? on I i co 

b89±158 


W1 o I 1 OO 


^70 O 1 1 O O 


n 1 c I i no 

915±198 


r*m I on a 

bUl±2U4 










U//1 / /Uo 


54bb3. 51324 


77.81 


46.50 


A A o _l_ 1 in 
442±11U 


a a o_l_ 1 1 o 
4b2±llo 


481±9U 


353±lb5 


5U3±119 


C7tr 1 ion 

b /5±139 


4Ub±19U 


959±2b5 


noi± A A 
922±41 


351±b81 


07/17/08 


54663.53845 


73.94 


48.31 


390±96 


417±96 


653±76 


479±106 


650±128 


511±147 


639±156 


944±294 


820±624 


322±998 


07/17/08 


54663.59463 


62.22 


49.86 


107±110 


281±98 


148±121 


117±172 


365±139 


673±158 


662±232 


902±261 


971±474 


913±793 


07/17/08 


54663.62064 


55.59 


48.99 


526±180 


250±151 


203±145 


318±174 


618±164 


584±219 


790±286 


1630±309 


1998±772 


1771±892 


07/18/08 


54664.50661 


78.33 


46.17 


498±150 


502±126 


683±120 


680±164 


756±152 


739±181 


946±266 


565±383 


1043±503 


1502±723 


07/18/08 


54664.54028 


73.14 


48.57 


750±176 


633±160 


626±157 


670±143 


494±170 


575±172 


846±212 


1027±518 


1245±591 


1506±790 


07/18/08 


54664.59729 


60.90 


49.79 


536±115 


527±85 


367±99 


495±115 


605±166 


878±140 


1049±171 


1179±241 


1788±447 


994±681 


07/18/08 


54664.61968 


55.12 


48.88 


188±123 


285±109 


230±116 


146±139 


505±125 


867±143 


390±212 


810±272 


1481±606 


2388±898 



Note. — Date format: month/day /year. MJD: Modified Julian Day = Julian Day - 2400000.5. B v : projected baseline in m. Az: baseline azimuth East of North in degrees. 
Ny. calibrated astrophysical null (x 10 5 ) at wavelength A. 
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Table 3: Summary of the measurements fitted in this study. The last column indicates which data are used 
for each run of simulations (A: ALL wavelengths from 2 to 13 ^m, S: SHORT wavelengths only, from 2 to 
11 /im, L: LONG wavelengths only, from 11 to 13 /im). [1] lLahuis et al.l (|2006l ) 



Wavelength (jum) 


Flux[Jy] 


Uncertainty (la) 


Instrument 


fit 


11.2328 


12.42 


0.27 


Spitzer/IRS (c2d) [1] 


— L 


11.7363 


11.37 


0.35 


Spitzer/IRS (c2d) [1] 


— L 


11.9366 


11.07 


0.26 


Spitzer/IRS (c2d) [1] 


-L 


Wavelength (/im) 


Null 


Uncertainty (lex) 


Instrument 


fit 


8.25 


0.008381 


0.003551 


KIN 2007 SB 


AS- 


8.74 


0.006274 


0.002116 


KIN 2007 SB 


AS- 


9.20 


0.007156 


0.002106 


KIN 2007 SB 


AS- 


9.75 


0.009150 


0.002196 


KIN 2007 SB 


AS- 


10.22 


0.005664 


0.002738 


KIN 2007 SB 


AS- 


10.68 


0.006289 


0.002175 


KIN 2007 SB 


AS- 


8.25 


0.006966 


0.002105 


KIN 2007 LB 


AS- 


8.74 


0.005859 


0.002104 


KIN 2007 LB 


AS- 


9.20 


0.008484 


0.002234 


KIN 2007 LB 


AS- 


9.75 


0.009504 


0.002671 


KIN 2007 LB 


AS- 


10.22 


0.008586 


0.002648 


KIN 2007 LB 


AS- 


10.68 


0.006852 


0.003473 


KIN 2007 LB 


AS- 


8.25 


0.003848 


0.001414 


KIN 2008 SB 


AS- 


8.74 


0.004776 


0.001414 


KIN 2008 SB 


AS- 


9.20 


0.003602 


0.001414 


KIN 2008 SB 


AS- 


9.75 


0.003854 


0.001414 


KIN 2008 SB 


AS- 


10.22 


0.005990 


0.001414 


KIN 2008 SB 


AS- 


10.68 


0.008624 


0.001414 


KIN 2008 SB 


AS- 


11.21 


0.008571 


0.002121 


KIN 2008 SB 


A-L 


11.72 


0.011752 


0.002475 


KIN 2008 SB 


A-L 


12.20 


0.014865 


0.003536 


KIN 2008 SB 


A-L 


12.69 


0.014790 


0.005657 


KIN 2008 SB 


A-L 


8.25 


0.008658 


0.001414 


KIN 2008 LB 


AS- 


8.74 


0.007888 


0.001414 


KIN 2008 LB 


AS- 


9.20 


0.008950 


0.001414 


KIN 2008 LB 


AS- 


9.75 


0.008292 


0.001414 


KIN 2008 LB 


AS- 


10.22 


0.008310 


0.001414 


KIN 2008 LB 


AS- 


10.68 


0.008367 


0.001414 


KIN 2008 LB 


AS- 


11.21 


0.009132 


0.002121 


KIN 2008 LB 


A-L 


11.72 


0.010642 


0.002475 


KIN 2008 LB 


A-L 


12.20 


0.011935 


0.003535 


KIN 2008 LB 


A-L 


12.69 


0.008260 


0.005657 


KIN 2008 LB 


A-L 


Wavelength (/im) 


Fractional Excess 


Uncertainty (la) 


Instrument 


fit 



2.18 



0.0088 



0.0012 



VLTI fractional excess AS- 



-24- 



Table 4: Parameter space explored with GRaTer 



Parameter 


Explored 


range 


Values 


Distribution 


r [AU] 


[0.05, .. 


.8.0] 


45 


log 


a 


[-9.0, . . . 


, 0.0] 


10 


linear 


flmin [yum] 


[0.01, ... 


, 100] 


15 


log 


H 


[-6.0, . . . 


-2.5] 


10 


linear 


V C 


[0, • • , 


19] 


20 


linear 


VSi 








10 




fixed 




a max [mm] 


1.0 




fixed 




M dust [M e ] 


> 




fitted 





Notes - We use the astron omical silicat e s from of iLi fc Greenbergl (119971 ) and the ACAR sample of 
carbonaceous material from IZubko et al.l (119961 ). Md us t is scaled independently when adjusting the 
surface density Eq at the peak position to fit the SED. 
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Table 5: Best-fitting parameters for the three approaches. The bracketed values correspond to the 
1-cr confidence intervals derived from our Bayesian analysis (possibly several families of solutions can 
co-exist), while the individual values are those corresponding to the smallest \ 2 on the grid. Some 
confidence intervals must be taken cautiously (semi-open intervals) because their margins correspond 
to the limits of the parameter space explored (see Fig. [5]). 



ALL SHORT LONG 



[9.6,19[ [8.3,19[ [0,19[ 

19 19 

[0.07, 0.14] U[0.33, 0.41] [0.08,0.11] [0.21, 0.62] U[0.88, 1.08] 

0.40 0.09 0.45 

]-9.0,-5.0] ]-9.0,-5.2] ]-9.0,-4.0] 

-9.0 -8.0 -9.0 

]0.01,0.08] ]0.01, 0.21] [1.5,69.0] 

0.02 0.01 6.6 

]-6.0,-5.3] ]-6.0,-5.2] ]-6.0,-3.6] 

-6.0 -6.0 -6.0 

[1.1,5.6] [1.5,2.4] [2.3,6.1xl0 5 ]U[l-7xl0 6 ,5.6xl0 6 ] 

2.2 2.2 75.5 



r (AU) 



a 



amm(/im) 



K 



M dust (10- 10 M e 



X 2 (dof) 



42.53 (28) 



29.04 (20) 



1.87 (5) 



